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Classical novae are potential sources of 7-rays, like the 
1.275 MeV gamma emission following 22 Na beta decay, that 
could be detected by appropriate instruments on board of fu- 
ture satellites like INTEGRAL. It has been shown that the 
production of 22 Na by novae is affected by the uncertainty 
on the 21 Na(p,7) 22 Mg rate and in particular by the unknown 
partial widths of the E x = 5.714, F = 2+, 22 Mg level. To 
reduce these uncertainties, we performed shell model calcu- 
lations with the OXBASH code, compared the results with 
available spectroscopic data and calculated the missing par- 
tial widths. Finally, we discuss the influence of these results 
on the 21 Na(p,7) 22 Mg reaction rate and 22 Na synthesis. 
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I. INTRODUCTION 

Classical novae are potential sources of 22 Na which 
beta decays (r = 3.75 years) towards the first excited 
state of 22 Ne followed by a prompt gamma ray emis- 
sion (£ 7 = 1.275 MeV) that could be observed 0. Even 
though it has not been detected yet § , this could happen 
with higher sensitivity instruments on board of future 
space missions like INTEGRAL (including a Ge based 
gamma-ray telescope ||) or gamma-ray focusing tele- 
scope projects j|. If conditions are favorable (i.e. an 
oxygen-neon nova at a distance of less than «2 kpc 
H) the 22 Na line could be detected § with the IN- 
TEGRAL spectrometer. Nova outbursts occur at the 
surface of an accreting white dwarf within a binary sys- 
tem. The accreted H-rich matter enriched with the C-0 
or O-Ne matter from the white dwarf undergo a ther- 
monuclear runaway that synthesizes new isotopes. The 
formation of 22 Na (in O-Ne novae) proceeds from ini- 
tial 20 Ne present in large quantities through the two 
possible paths : 20 Ne(p,7) 21 Na(/3 + ) 21 Ne(p,7) 22 Na and 
20 Ne(p,7) 21 Na(p,7) 22 Mg(/3+) 22 Na. The first path has 
been found JtJ to be more favorable to 22 Na formation 
because of its longer time scale. The preferred path 
is governed by the competition between the 21 Na f3 + ~ 
decay and the 21 Na(p,7) 22 Mg reaction whose rate re- 
mains uncertain mainly because of the unknown reso- 
nance strength associated with the E x = 5.714, J T = 
2+ 22 Mg level §. 

Estimates of the 21 Na(p,7) 22 Mg reaction rate @] have 
been provided |j) 10|, considering the first three levels lllll 



above the proton threshold (5.501 MeV). With the excep- 
tion of the first level, the total widths can be identified 
with the proton widths, so that W7 rj ujTj. The reso- 
nance strength associated with the first E x =<3.714 MeV, 
J 7r =2 + level suffers from a significant uncertainty that af- 
fects the 21 Na(p,7) 22 Mg rate in the temperature domain 
of nova nucleosynthesis 0. For this level, only the total 
width is known experimentally (r=16.5±4.4 meV [|Tl")) 
and to calculate the corresponding resonance strength, 
estimated proton widths |9|,|l0| have been used together 
with the relation r = T p + T 7 . The two estimates (l p = 
and 6> 2 = 0.01 fj or l p = 2 and 6» 2 = 0.5 fuj) lead to simi- 
lar values (7= 3.4 or 3.8 meV) very close to the maximum 
value (r/4), obtained when r 7 = T p = Y/2. However, 
based on the data available for the 22 Ne £^=6.120 MeV 
analog level jL2| , it was argued H that the proton width 
could be much smaller because i) the radiative width es- 
timated from the analog level is such that r 7 w V and 
ii) the neutron spectroscopic factor in the analog level 
should be very small according to experimental data [ fl3| . 
Accordingly, values of ui-f = 2.5, 0.25, 0.0 meV have been 



adopted for upper (F 1 



r/2), recommended (with 



the usual 0.1 reduction factor) and lower limit for the 
calculation of 22 Na production in novae 0. This induces 
a factor of 10 uncertainty on the rate around a temper- 
ature of » 10 s K, typical of novae, and a factor of up 
to 3 in the 22 Na yields @. Another, much less impor- 
tant, source of uncertainty |7| comes from the assumed 
value ||[l0) for the radiative width of the third (E x ; J 71 
= 5.965 MeV; + ) level. It is important to reduce this 
uncertainty on the 22 Na yield that directly affects the 
detectability distance of the 1.275 MeV gamma emission 
in order to interpret future nova observations. In conse- 
quence, we performed shell model calculations of spectro- 
scopic factors and radiative strengths for 22 Ne and 22 Mg 
nuclei. In this paper, we first compare calculated values 
with existing experimental spectroscopic data in order to 
validate the calculations. We re-analyze existing exper- 
imental data to extract supplementary information on 
missing spectroscopic factors. From this analysis, we de- 
rive better estimates for the spectroscopic factor of the 
5.714; 2+ level and the radiative width of the 5.965 MeV; 
+ level. Finally, we discuss the influence of these new 
values on the 21 Na(p,7) 22 Mg rate and 22 Na production 
in novae. 



{E x \ ,r = 5.714; 2+, 5.837; (0-5) and 5.965 MeV; 0+) 
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II. SHELL MODEL CALCULATIONS 

In order to estimate spectroscopic factors and radiative 
widths, we have performed shell model calculations using 
the OXBASH code [[uj. Since we are interested only in 
positive parity states, we used the well-known USD in- 
teraction of Wildenthal for the sd shell model space. 
The results of the calculations in comparison with the ex- 
perimental spectra of 22 Ne and 22 Mg are shown in Figs. 0, 
H and Table |. The correspondence of the experimental 
and theoretical results is remarkably good. 
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from the scattering phase shifts in the Woods-Saxon po- 
tential with the depth required to reproduce the experi- 
mentally known energy of the resonance. The contribu- 
tion of / = 2 transfer to the 5.714 MeV state in 22 Mg is 
negligible as compared to I = transfer and we obtain 
the values of T 0d5/2 = 2 x 1(T 6 eV, T Qda/2 = 2x 10" 7 eV, 
T 0s = 4.5 x 1CT 3 cV, that leads to T p = 4.5 meV. 
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FIG. 1. Positive parity energy levels of T=l A=22 nuclei 
calculated with the USD interaction in comparison with the 
experimental level schemes of 22 Ne and 22 Mg up to 5 MeV. 

In Table | the calculated neutron spectroscopic fac- 
tors in 22 Ne are compared with the experimental val- 
ues obtained Jl3[ through the neutron stripping reaction 
21 Ne(d,p) 22 Ne. The agreement between calculated and 
experimental values is very good except for the levels la- 
beled (in ref. Jll|) 6.120, 6.350; 6 + (no experimentally de- 
termined spectroscopic factors) and 7.341; + . However, 
these discrepancies can be explained when new spectro- 
scopic data |l2| are included (see the following section). 

Assuming the equality between spectroscopic factor of 
conjugate reactions and according to the good agreement 
between our calculations and experiment we can confi- 
dently use them to obtain the 22 Mg proton width of the 
E x ; J* = 5.714; 2+ level. The fourth calculated 2+ 
state at 6.179 MeV corresponds to the fourth 2 + state 
of 22 Ne at 6.120 MeV and to the 2+ state of 22 Mg at 
5.714 MeV (Fig. ||) which are of main interest here. To 
check the correctness of the assignment and the proxim- 
ity of the calculated and physical 2+ state, we compared 
calculated radiative strength with those available experi- 
mentally Jl6| for the £^=6.120 MeV level in 22 Ne. There 
is a fair agreement between the experimental and theoret- 
ical values as shown in Table || increasing the confidence 
in the assignment. The calculated spectroscopic factors 
for the 2 + state are small as expected from experiments 
on 22 Ne (see Table Q). The corresponding proton reduced 
widths are obtained, using the relation T p = C 2 S T s , p , 
where the single-particle width r s p . has been estimated 
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FIG. 2. Same as Fig. |l| but in the energy range between 5 
and 6.8 MeV. In addition, for 22 Mg, the threshold for proton 
emission is shown (dashed line) together with the position of 
the Gamow peak for temperatures relevant to nova outbursts 
(arrows). 

We also obtained the radiative width of the third level 
above threshold (E x ; J 7 ' = 5.965 MeV; 0+) by calcu- 
lating B(E2) or B{M1) for the transition to the lower 
lying 2 + and 1 + levels (Table [Hi]). The calculated value, 
T 7 = 33.6 meV leads to a resonance strength of u>~f = 
4.2 meV, only slightly higher than the estimated value 
|i of 2.5 meV. 
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III. RE ANALYSIS OF EXPERIMENTAL DATA 

We performed a DWBA analysis of some of Neogy et 
al. Jl|] data (6.120, 6.350 and 7.341 levels) to extend the 
comparison between calculated and experimental values 
and to put more constraint on the 6.120 MeV 22 Ne (or 
5.714 MeV 22 Mg) level spectroscopic factor. For this pur- 
pose, we used the ECIS code jlTj with the same optical 
potential parameters as Neogy et al. 

The 6.350 MeV level (Fig. 5 of Neogy et al. fn|) was 
assumed to have a 6 + spin and parity and accordingly 
(I = 4) no DWBA analysis was performed at that time. 
However, since then [Q, a 6 + level has been located at 
6.311 MeV and a 4+ level at 6.345 MeV. Hence, it is 
almost certain that the 6.350 MeV level in ref. jl3| was 
unduly identified with the 6 + instead of the 4^1cvcl. 
Accordingly, we made a DWBA analysis of the Neogy et 
al. data assuming that the reported level is a 4 + and 
extracted a spectroscopic factor. This new value is close 
(within a factor of two) to the calculated, 4^" level, value 
(see Table |). Hence, shell model calculations reinforce 
the idea that the 6.35 MeV level seen by Neogy et al. JljJ 
is the 6.345; 4 + one instead of the 6.311; 6 + . 

Our shell model calculations lead to a very small spec- 
troscopic factor for the 7.341; + level in complete con- 
tradiction with the value reported by Neogy et al. How- 
ever, less than three keV above (E x — 7.344 MeV) lies a 
J v = (3,4) + level whose calculated spectroscopic factors 
(for the 3;j~ and 4^ states) agree much better with those 
extracted from experimental data (Table |). Hence, the 
experimentally determined spectroscopic factor could be 
attributed to the 7.344 MeV; (3,4)+ level rather than to 
the 7.341; 0+ one. 

In order to put more constraint on the 6.120 MeV 22 Ne 
level we also performed a DWBA analysis of the Ne- 
ogy et al. data for this level. (This analysis was not 
performed in the original work |13j| because "the angu- 
lar distribution does not exhibit characteristics of direct 
reactions".) The experimental data and DWBA cross 
sections are represented in Fig. ||. As expected, the cal- 
culated transfer cross sections are more forward peaked 
than the experimental angular distribution suggesting a 
strong contribution from fusion reactions. When the the- 
oretical spectroscopic factors are used, the d 5 / 2 contribu- 
tion is negligible while the Si/ 2 contribution is compatible 
with experimental data except for the most forward an- 
gle. Requesting that the DWBA cross sections remain 
below all experimental data points lead to upper limits 
for the spectroscopic factors of C 2 S ^ 0.0025 or 0.015 
for / = or 2 respectively. While the I = 2 upper limit is 
fully compatible with shell model calculations, the I = 
one is a factor of four below the calculated spectroscopic 
factor. From this upper limit, we obtain T p 1 meV. 
One can note however that this limit should be taken 
with caution as it comes from a single data point at the 
smallest angle. As it can be seen in Fig. 1 of Neogy et al., 
the 6.120 MeV, 22 Ne peak is close to an other one from 



23 Ne (arising from a (d,p) reaction on 22 Ne in the tar- 
get.) Resolving these two peaks at a lower angle should 
be more difficult because of the unfavorable evolutions of 
both their energy separation and relative heights. Hence 
one cannot exclude that the experimental error bars were 
underestimated in this case. 
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FIG. 3. Angular distribution for the 21 Ne(d,p) 22 Ne neu- 
tron transfer reaction on the 6.120 MeV; 2 + level. The ex- 
perimental data is from Neogy et al. |l2} ]. The curves repre- 
sent DWBA calculations with calculated spectroscopic factors 
(dashed lines) or adjusted spectroscopic factors (solid line). 



IV. IMPLICATIONS ON THE 21 Na(p, 7 ) 22 Mg RATE 

The contribution of the E x ; J 71 — 5.714; 2+ level to the 
21 Na(p,7) 22 Mg rate depends directly from its adopted 
proton width. From experimental data ]l3[ ] we deduced 
an approximate upper limit (see comment above) for this 
width (r p 5s 1 meV). Our shell model calculations give a 
slightly higher value (r p = 4.5 meV). Using the experi- 
mental total width (r = 16.5 meV), one obtains 7 « 3.2 
or 1 meV, and lj 7 — 2. or 0.6 meV, when using the calcu- 
lated or experimental upper limit for the proton widths. 
These values are close to the first (3.4 @ and 3.8 meV 
|]l0| ) and more recent estimated strengths (2.5 upper and 
0.25 meV nominal values [Q). One important conclu- 
sion resulting from the shell model calculations is that 
it is most unlikely that the spectroscopic factor is much 
smaller than one meV. This would exclude the lower limit 
for the rate |?J obtained with a null spectroscopic factor. 

The contribution of the 5.965 MeV; + level is obtained 
through our shell model calculation of its radiative width 
= 33.6 meV and w 7 = 4.2 meV). As in previous works 
(H we assume that the 5.837 MeV; (0-5) level, is the 
analog of the 5.910; 3~ level in 22 Ne. This assignment 
made by flllfl , is not present anymore in |12f] but is likely 
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from the examination of Fig . || and because of their simi- 
lar gamma decay modes |16|| . Accordingly, we take T 7 = 
13 meV from the analog 22 Ne level and oj-f = 11.4 meV. 
The direct capture term ||,|| is also left unchanged as it 
is based on experimental spectroscopic factors. 

Following the good general agreement between cal- 
culated and experimental quantities, we derive the 
21 Na(p,7) 22 Mg rate using the shell model calculated val- 
ues presented above. The resonant part of the rate is 
approximated as usual by J^i^i exp(— Sj/Tg) with A; 
= 334., 1862., 686. for B, = 2.52, 3.95, 5.49 respec- 
tively. The resulting rate is very close to the previous 
rates P,^0|,p[ and of the recent upper rate limit j?j] but 
it is now put on a safer ground, in the domain of nova 
nucleosynthesis, as it now relies on shell model calcu- 
lations rather than estimates. Using the experimental 
upper limit instead of the calculated one for the first res- 
onance strength would lead to a rate lying within a factor 

and from the nom- 
(The proximity 
of these different rates is due to the strong experimental 



of three of the rates from refs. ||9JJ1 
inal and upper rate limit from ref. 



r, 



16.5±4.4 meV 



21 Na(p, 7 ) 22 Mg 
yields in no- 



constraint introduced by T 
0") 

In conlusion, the lower limit for the 
rate, used for calculating JtJ 22 Na 
vae seems now excluded. Unfortunately it was also 
the more efficient for 22 Na production through the 
21 Na(/?+) 21 Ne(p,7) 22 Na chain. Hence the highest 22 Na 
yields reported Q are now excluded. A precise conclusion 
on gamma emission detectability will require further hy- 
drodynamical calculations of nova outbursts. However, 
the 21 Na(p,7) 22 Mg rate resulting from this analysis is 
not too far from the nominal rate used in previous cal- 
culations j?J so that the nominal detectability distance 
of 22 Na gamma emission || should not be too much af- 
fected. 
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TABLE I. Experimental, positive parity, energy levels of Mg and Ne and spectroscopic 
factors deduced from the 21 Ne(d,p) 22 Ne reaction compared with calculated (T = 1 A — 22) values. 
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a E x and J 71 from Ref. Q 
b E x and J" from Ref. H 

c (2 J + 1)5 from Ref. ||lo| unless othewise stated. 
d Our analysis of Neogy et al. data 0|. 

e Assuming that the 6.35; 6 + 22 Ne data in Neogy et al. Jl3| corresponds to the 6.345; 4 + 
level. 

Assuming that the 7.341; + 22 Ne data in Neogy et al. (ll| corresponds to the 7.344; 
(3,4)+ level. 



TABLE II. Experimental and theoretical reduced transi- 
tion probabilites from 6.120 MeV 2+ state 22 Ne. 

Transition Experimental a Theoretical 

2+ -> 0+ B(E2) = (1.91 ± 1.11) e 2 fm 4 B{E2) = 3.1 e 2 fm 4 

2+ 2+ B(E2) = (1.15 ± 1.15) e 2 fm 4 B(E2) = 1.8 e 2 fm 4 

2+ 2+ B(M1) = (0.046 ± 0.023) fj.% B(M1) = 0.138 fj.% 

a From Ref. @ 
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TABLE III. Theoretical reduced transition probabilites 
from Oj state of which is assumed to correspond to the 5.965 
MeV 0+ level in 22 Mg. 

Transition Transition rates 

Oj -> 2+ B(E2) = 16.92 e 2 fm 4 

0+ -> 2+ B(£2) = 8.0 e 2 fm 4 

0+ 2+ B(£2) = 6.2 e 2 fm 4 

0+ -> 2+ B(E2) = 0.01 e 2 fm 4 

Oj -> 1+ B(M1) = 0.57 pif 
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